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We show that normal fluid eddies in turbulent helium II polarize the tangle of quantized vortex 
lines present in the flow, thus inducing superfiuid vorticity patterns similar to the driving normal 
fluid eddies. We also show that the polarization is effective over the entire inertial range. The results 
help explain the surprising analogies between classical and superfiuid turbulence which have been 
observed recently. 
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Our concern is the experimental evidence that turbulent helium II appears similar to classical hydrodynamic 
turbulence For example, Smith et al || found that the temporal decay of helium II turbulence behind a towed 
grid is the same as that expected in an ordinary fluid. Maurer & Tabeling |3| observed the classical Kolmogorov k~ 5 / 3 
dependence of the energy spectrum on the wavenumber k in helium II continuously agitated by rotating propellers. In 
other experiments it was found that if helium II is forced at high velocity along pipes and channels the same pressure 
drops Q are detected which are observed in a classical liquid. Furthermore, when a sphere moves at high velocity in 
helium II, the same drag crisis || is measured that occurs if the fluid is air. 

These apparently classical results are surprising because helium II is a quantum fluid. According to Landau's 
two-fluid theory, it consists of the intimate mixture of an inviscid superfiuid component and a viscous normal fluid 
component. The latter is similar to a classical Navier - Stokes fluid, so, when made turbulent, it consists of eddies of 
various sizes and strengths. On the contrary quantum mechanics constrains the rotational flow of the superfiuid to 
discrete quantized vortex lines, each vortex characterized by the same quantum of circulation T = 9.97 x 10~ 4 cm 2 /sec. 
Unlike what happens in a classical Euler fluid, superfiuid vortex lines can reconnect with each other 0. They also 
interact with the normal fluid via a linear mutual friction force |Q, and, when helium II is made turbulent, they form 
a disordered, apparently random tangle. Clearly the experimental results described above call for an explanation 
in terms of the basic physical ingredients of the problem (normal fluid eddies and superfiuid vortices) and their 
interaction. What is remarkable is that these classical aspects of helium II turbulence are observed to be independent 
of temperature, whereas the relative proportion of superfiuid and normal fluid is a strong function of temperature. 

A possible solution of the puzzle is that both fluids independently obey the classical Kolmogorov scaling. A few 
years ago Nore et al || solved the Gross-Pitaevskii equation for a tangled system of vortices and noticed that the 
energy spectrum follows a fc -5 / 3 dependence. Unfortunately the effect only took place over a small range in k space 
and further calculations performed using the vortex filament model Q did not reproduce it. Another difficulty with 
this argument is that one has also to explain why the mutual friction coupling between the two fluids does not change 
the shape of the spectra which the two fluids would have independently of each other. 

A second explanation of the experiments has been suggested by Vinen jl(| . He argued that at spatial scales larger 
than the average separation S between the quantized vortex lines, the superfiuid and normal fluid are coupled by 
a small degree of polarization of the almost random tangle of superfiuid vortex lines. If that is the case, on these 
scales helium II behaves as a single fluid of density p = p s + p n , consistently with the experiments (p s and p n are the 
superfiuid and normal fluid density respectively). 

The aim of this letter is to support Vinen's argument with quantitative evidence of polarization. First we shall 
introduce some simple models which, although very idealized, capture the essential physics of polarization. Secondly 
we shall look for evidence of polarization by direct numerical simulation. 

Our first model is concerned with the reaction of superfiuid vortices to a normal fluid shear. Consider a row of 
point vortices of alternating circulation ±r initially set along the x axis at distance 5 from each other. We assume 
that the normal fluid is v„ = V n cos(fcy)x. The governing equation of motion of a vortex point is [ fl3f 

dy 

-- = ±aV n cos(ky), (1) 
at 

where a is a known temperature dependent friction coefficient . The solution of Eq. ([!]) corresponding to the initial 
condition y(0) = is 

^ ) = l(~I + ton " 1(e±QW " t) )- (2) 

Given enough time (t — > oo), positive and negative vortices will reach stable locations yoo = ±ir/2k respectively. In 
a turbulent normal flow, however, the shear does not last longer than few times the turnover time r ~ l/oJ n ~ l/kV n . 
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Since a is small (it ranges from 0.037 at T = 1.3K to 0.35 at T = 2.16K) we have y(r) « ±6 where b = a/k. Within 
the lifetime of the shear we have thus created a separation 2b between positive and negative vortices, that is to say we 
have polarized the initial configuration. The velocity of this Karman - vortex street is approximately V s « ttT/2S 
for kS << 1 in the direction along the x axis where the normal fluid (which induced the polarization in the first 
place) is stronger. The result suggests that it is not necessary to create extra vortex lines to generate a superfluid 
pattern that mimics the normal fluid one - rearranging existing vortices is enough. We also notice that the induced 
polarization is proportional to a. 

Our second model is concerned with the expansion of favourably oriented superfluid vorticity. We think of the 
superfluid vortex tangle as a collection of vortex rings of radius approximately equal to the average separation of 
vortices in the tangle, Rq ~ <5. We assume for simplicity that the rings are on the x,y plane with equal numbers of 
rings oriented in the ±z directions. Depending on whether they have positive or negative orientation, the rings move 
along ±z with self induced speed given by 

TC 

where C = In (8i?o/ao) — 1/2 is a slowly varying term and a w 10~ 8 cm is the vortex core radius. Now we apply a 
normal fluid velocity V n in the z direction. The radius R of a ring is determined by 

^ = ^(V n -Vn), (4) 
at p s L 

where 7 is a known friction coefficient and 7/p s T « a at almost all temperatures of interest. Eq. (|]) shows that V n 
selectively changes radius and velocity of vortex rings moving in opposite directions. A ring which grows (shrinks) 
by an amount SR = aSt(V n — Vr ) in time St slows down (speeds up) by an amount 5Vr = Vr 5R/ Rq. In this way 
a superflow is induced in the same direction of the normal fluid which induced the polarization in the first place. A 
simple estimate of the spatial averaged magnitude of this superflow yields p| V s w 3Vr 5R/ Rq. 

Our third model is concerned with the rotation of existing superfluid vorticity, We represent a superfluid vortex 
line as a straight segment pointing away from the origin and study how its orientation is changed by a normal fluid 
rotation about the z axis. Using spherical coordinates (r,9,(j>), we assume that the vortex is initially in the plane 
8 = 7r/2. The normal fluid's velocity is v n = (0, 0, QrsinO) and the motion of the vortex segment is determined by 
@ 

d9 

— = -an sin{6), (5) 

and dr/dt = and d<j>/dt — 0. The solution to Eq.(||) is 

6{t) = 2 tan-\e- am ), (6) 

with r and 4> constant. Given enough time, the vortex segment will align along the direction of the normal fluid 
rotation [6 — ► for t — > 00), but the lifetime r of the eddy is only of the order of r « 1 /O, so the vortex can only 
turn to the angle 

0(t) w tt/2 - a. (7) 

Despite the smallness of the angle, the effect is sufficient to create a net polarization of the tangle in the direction 
of the normal fluid's rotation, provided that there are enough vortices. The following argument shows how this is 
possible. The normal fluid is like a classical viscous Navier - Stokes fluid, and, if left to itself, its spectrum would 
obey Kolmogorov's law 

E k = Ce 2/3 k~ 5/3 , (8) 

Eq. (||) is valid in the inertial range l/£o < k < I/77 in which big eddies break up into smaller eddies, transferring 
energy to higher and higher wavenumbers without viscosity playing a role. Here k is the magnitude of the three 
dimensional wavevector, e is the rate of energy dissipation per unit mass, C is a constant of order unity, £q is the 
integral length scale (the scale at which energy is fed into the energy cascade) and ij is the Kolmogorov scale at which 
kinetic energy is dissipated by the action of viscosity. In reality the normal fluid is not alone but is forced by the 
quantized vortex filaments. We know little of the effects of this forcing (it has been studied only for very simple 
geometries [^5|) so, for lack of further information, we assume that the classical relation Eq. (||) is valid for the normal 
fluid. 
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The quantity which is often used to describe the intensity of the superfiuid vortex tangle is the vortex line density 
L (length A of vortex line per volume V) because it is easily measured by detecting the attenuation of second sound. 
From L one infers the average separation between vortices, 5 « L~ x / 2 . The quantity TL can be interpreted as the 
total rms vorticity of the superfiuid. Note that the net amount of superfiuid vorticity in a particular direction can be 
much less than TL (even zero, if the tangle is randomly oriented). 

The key question is whether, as a result of mutual friction, sufficient quantized vortex lines can re-orient themselves 
within a normal fluid eddy of wavenumber k so that the resulting net superfiuid vorticity matches the vorticity ujk of 
that eddy. The process must take place in a time scale shorter than the typical lifetime of the eddy, which is of the 
order of few times the turnover time l/u>k- At this point we use the result of the third model, for which the initial 
condition 9(0) = tt/2 represents the average case. In doing so we remark that, if the initial orientation of the vortex 
is toward the origin rather than away from it, then the vortex segment turns to 9(t) rs 7r/2 + a rather than tt/2 — a 
but still contributes to positive vorticity in the z direction. Therefore in the time re-ordering of existing vortex 

lines creates a net superfiuid vorticity oj s of the order of aLT /3 in the direction of the vorticity Uk of the driving 
normal fluid eddy of wavenumber k. Since Uk is approximately LUk = \/WE\ we have uj k = C 1/2 e 1/3 fc 2/3 . Matching 
uj s and Wk would then require 

icTT > C^V^k 2 / 3 . (9) 

The normal fluid vorticity increases with k and is concentrated at the smallest scale (k f=a 1 / 77) , so a vortex tangle 
with a given value of L may satisfy the above equation only up to a certain critical wavenumber k c . Substituting 
e = v^/rf where v n is the normal fluid's kinematic viscosity (the viscosity of helium II divided by p„), we obtain 
(S/ri) = C- 1 / 4 (a/3) 1 / 2 (r/i/„) 1 / 2 (? 7 fe c )- 1 / 3 . If k c « I/77 then 



In the temperature range of experimental interest T/v n ranges from 0.43 at T — 1.3K to 5.86 at T = 2.15K, so 
S/rj = 0(1) and we conclude that matching of the two vorticities (k c « 1 / 77) is possible throughout the inertial range. 

Due to the computational cost, numerical simulations of superfiuid turbulence do not produce vortex tangles dense 
enough to cover the range k < l/£ and determine unambiguously the dependence of the energy spectrum on k 
in this range. To make progress in the problem and confirm the above arguments we study the reaction of the 
superfiuid vortex tangle to a single scale ABC normal flow jl6| given by v n = (A sin (kz) + C cos (ky), B sin (kx) + 
A cos(kz),C sm(ky) + B cos(kx)) where k = 2tt/X is the wavenumber, A is the wavelength and A, B and C are 
parameters. ABC flows are solutions of the Euler equation and the forced Navier - Stokes equation and have been 
used as idealized model of eddies in fluid dynamics, magneto-hydrodinamics and superfiuid hydrodynamics [ |L7| . For 
the sake of simplicity we take A = B = C and A = 1. 

Following Schwarz jlEfl , we represent a superfiuid vortex filament as a space curve s = s(£,i) where £ is arclength 
and t is time. Neglecting a small transverse friction coefficient, the curve moves with velocity 

ds 

— = v si + as' x (v„ - Vn). (11) 

where s' = ds/d£ and the self induced velocity v si is given by the Biot - Savart integral 

T f (r — s) x dr 
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The calculation is performed in a cubic box of volume V — 1 cm 3 with periodic boundary conditions. The numerical 
technique is standar d |lq] and the details of our algorithm, including how to perform vortex reconnections, have been 
published elsewhere |17| . 

We start the calculation with N — 50 superfiuid vortex rings set at random positions and orientation and integrate 
in time at a variety of temperatures (a = 0.1, 0.5 and 1.0) and normal fluid's velocities (A = 0.01, 0.1, 1.0 and 
10.0 cm /sec). The vortex length (A = 76.8 cm at t — 0) increases or decreases depending on whether the ABC flow 
is strong enough to feed energy into the normal fluid via instabilities of vortex waves (for example, for a = 1.0, the 
final length A is as high as 781.3 cm at A = 10.0 cm/sec, and as low as 56.67 cm at A = 0.01 cm/sec) The rings 
interact with each other and with the normal fluid, get distorted, reconnect, and soon an apparently random tangle 
is formed (see figure 1). 

The quantity < cos (0) >=< s' • u) n >, which we monitor during the evolution, gives us the tangle - averaged 
projection of the local tangent to a vortex in the direction of the local normal fluid vorticity, <l) n = (l/cj n )a; n , where 
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u n — V x v n . At t = < cos (8) >= due to the random nature of the initial state, and it is apparent from figure 2 
that < cos (9) > increases with time, no matter whether A decreases or increases. 

The results are analyzed in figure 3. From the simple models descrived above we expect that the polarization 
induced by the normal fluid vorticity is proportional to a. We also know from the discussion above that we should 
restrict the analysis to times t < t where r = l/u>„ with ui n = ^/3Ak is the lifetime of the normal fluid eddy, which 
we assume to be the same as the eddy's turnover time. It is apparent from the figure that, no matter whether the 
tangle grows or decays, approximately the same polarization takes place for t/r < 1. 

In conclusion we have put the theory of superfluid turbulence on firmer ground. Using simple models which capture 
the essential physical mechanisms of polarizaton and then a numerical simulation, we have shown that, within the 
lifetime of a normal fluid eddy of wavenumber k, superfluid vortex lines can rearrange themselves so that the superfluid 
vorticity and the normal fluid vorticity arc aligned. Provided that enough vortex lines are present, vorticity matching 
should take place over the entire inertial range, up to wavenumbers k of the order of l/£. 

Our result has theoretical and experimental implications. Numerical simulations of vortex lines driven by normal 
fluid turbulence |ll| show a superfluid energy spectrum in the accessible region k > l/S. More intense (hence com- 
putationally expensive) vortex tangles should be investigated to explore the region k << l/S where, as a consequence 
of our result, we predict the classical fc~ 5 / 3 dependence. Clearly an important issue which must be investigated is the 
nonlinear saturation of the polarization process. On the experimental side, our result supports the use of helium II 
to study classical turbulence. This has been done recently by Skrbek et al p3| who exploited the physical properties 
of liquid helium to study the decay of vorticity on an unprecedented wide range of scales. 
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FIGURE CAPTIONS 
Figure 1. Vortex configuration at t — 0.123 sec for a = 0.5 and A = 1.0 cm/ sec. 
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Figure 2. Average polarization < cos (8) > versus time t computed for different values of A and a. 

Figure 3. < cos (9) > /a versus scaled time t/r where r is the eddy's lifetime (same symbols as in Figure 1). 
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